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We discretize the Hamiltonian scalar constraint of three-dimensional Riemannian gravity on 
a graph of the loop quantum gravity phase space. This Hamiltonian has a clear interpretation 
in terms of discrete geometries: it computes the extrinsic curvature from dihedral angles. 
The Wheeler-DeWitt equation takes the form of difference equations, which are actually 
recursion relations satisfied by Wigner symbols. On the boundary of a tetrahedron, the 
Hamiltonian generates the exact recursion relation on the 6j-symbol which comes from the 
. Biedenharn-Elliott (pentagon) identity. This fills the gap between the canonical quantization 

^j^' and the symmetries of the Ponzano-Regge state-sum model for 3d gravity. 

<_ 

I. INTRODUCTION 

Q ' Three-dimensional Riemannian gravity is now a quite well understood theory, much easier than the four- 
O^. dimensional case since it does not carry local degrees of freedom. In spite of its apparent simplicity, it 
5—1 has been a good testing ground to build new methods aiming at a background independent quantization 
of gravity. In particular, the quantization programs of loop quantum gravity and spin foams have been 
successfully applied, on the kinematical aspects of quantum geometry [H,!!] as well as on its dynamics |3|-t9i]. 

Still, there remain two important questions in 3d gravity, one in the canonical loop quantum gravity 
(LQG) approach, and the second relating to the spin foam formalism. Here, we will propose a common 
solution to both, showing they are two aspects of the same problem. To say things briefly, a direct and natural 
quantization of the scalar constraint in 3d gravity, which would reproduce known results, is still missing (see 
• [3] though, and the discussion below), while in the Ponzano-Regge spin foam approach, recurrence relations 
on the basic amplitudes, coming from the topological symmetry of the model, still need to be generated 
. from the canonical formalism. 
. . ' Through the paper, we will consider pure three-dimensional gravity as a Chern-Simons-BF theory. This 
. ^ ' is a topological gauge theory whose reduced phase space is the cotangent bundle to the moduli space of flat 
^ G connections (here with G = SU(2)). Such theories (like two-dimensional Yang-Mills, three-dimensional 
Chern-Simons, and BF in higher dimensions) display some intriguing relations with the representation theory 
of the local symmetry group when one wants to quantize before imposing the symmetries. Typically, one 
can write their partition functions as state-sums based on the 6j-symbol (or higher Wigner symbols, or 
g-deformed versions). Then the symmetry which makes the model topological is described by algebraic 
relations on the 6j-symbols, such as the Biedenharn-Elliott (BE), pentagon identity. 

Usually, the canonical quantization in the Hamiltonian framework is a preferred setting to study the 
symmetries. Hence one should be able to derive such algebraic characterizations from the generators of the 
symmetries in the Hamiltonian formalism. As for three dimensional gravity itself, the framework of LQG 
seems very suitable since the boundary data of the corresponding state-sum model (the Ponzano-Regge 
model [4I]) are the kinematical quantum numbers of LQG (on a graph dual to a triangulation) . This is the 
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reason why a sp ecific project of studying tlie symmetries of simple models from the LQG point of view was 



initiated in [111 ] using recursions on recoupling Wigner symbols and generalizations. The present paper is a 
natural follow-up which goes further in the 3d gravity case. Indeed, in spite of the well-established fact that 
the moduli space of flat SU(2) connections can be reached from SU(2) recoupling theory [3, @, @, 14] (and 



more recently [13])) there has been no derivation of the key recursion on the 6j-symbol from Hamiltonian 
generators in LQG (but see an interpretation in [l^). 

A.iU) l^'^' f H + AM) If f fl + ^-iO) f f| = 0. (1) 



This is a second order difference equation (whose coefficients will be given in the main text) which can be 
derived from the BE identity and in turn fully determines the 6j-symbol (lol. E3|. 

The above discussion should be put in parallel with the fact that we do not know the Wheeler-DeWitt 
equation in 3d LQG, 

H\i^)=0, (2) 

i.e. the quantum equation for the Hamiltonian scalar constraint H in the spin network basis. Spin network 
states in LQG are labeled by spins attached to links of embedded graphs. Hence, if the Hamiltonian is 
not graph changing (as expected in three dimensional gravity), the Wheeler-DeWitt should be a difference 
equation, with non-trivial coefficients. It should be possible to solve it so as to reproduce results known 
from other approaches, and make contact with the symmetries encoded in spin foams. The recursion ([T|) is 
obviously the natural candidate, when the 6j-symbol is the physical state, like on the 2-sphere triangulated 
with four triangles (details in the main text). 

We consider in the present note a classical quantity H^j on the phase space of LQG dual to a triangulation 
and turn it into an operator H^j which enjoys the following nice properties. 

• Hjjj can be understood as a discretization of the scalar Hamiltonian constraint and forms a first class 
system (we only checked that on a specific triangulation of the 2-sphere) . 

• It has a nice, natural interpretation in terms of discrete geometries and dihedral angles of flat simplices. 

• Hyj acts on spin network functionals, and is labeled by a vertex u of a spin network graph and a 
cycle containing it - this cycle bounds the face /. 

The main result is that Hyj is the missing link between the scalar constraint in LQG and the BE identity. 



The equation H^j {ip) = 0, when written in the spin network basis, produces difference equations of 
order 2. 

For trivalent faces, the difference equation is ([T|) and hence solved uniquely by 6j-symbols, which 
proiect the curvature around the face to zero, as desired (for the solutions which we expect more 
generally, see fiQ, and correspondingly some expected recursions can be found in 0,0). 

Even if we cannot fully solve the model with our new proposal for the Hamiltonian yet (because the 
case of faces of more than three links requires recursion^elations on arbitrary Wigner symbols), it is worth 
comparing with the previous LQG quantization. In [6!], the authors use a simpler form of the flatness 
constraint which makes the quantization easier, but is far from being applicable to four dimensional gravity. 
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Nevertheless, the contact with SU(2) recouphng is natural is this case (see also [Sj] which is very close), 
though not complete with respect to our proposal. A completely different proposal was made in to 
quantize directly the scalar constraint. The author made use of what is now known as the Thiemann's trick, 
whose advantage is that it can be straightforwardly extended to four dimensions, and a major success is 
certainly that it leads to an anomaly-free algebra. However diffeomorphisms of the canonical surface and the 
scalar constraint are not treated on equal footing, in contrast with what happens in state-sum models. We 
will argue that our discretization of the scalar constraint actually and surprisingly contain all the constraints 
on the graph. This is certainly the deep reason why the quantization of [il does not exhibit any clear link 
with the recoupling theory and the BE identity. 

In section [nl we review some elementary facts about 3d Riemannian gravity, its loop quantization, and 
the physical state on the 2-sphere. Then, in section IIIIl we present our proposal for the Hamiltonian 
operator on spin networks, describe its classical geometric interpretation and show it is first class (on a 
specific triangulation) . In section \W\ we proceed to the quantization, and prove the relation to the flatness 
constraint on trivalent faces. Finally, in section IVj we discuss higher order difference equations, coming from 
lifting our Hamiltonian to representations of higher spins. 



II. REVIEW OF 3D LOOP QUANTUM GRAVITY 
A. Canonical analysis of 3d gravity 

Three-dimensional Riemannian gravity can be formulated with a connection A on a principal SU(2)- 
bundle P over a three-manifold M, and a cotriad e, that is an Ad(P)-valued 1-form. Locally, the connection 
is seen as a Lie algebra su(2)-valued 1-form over M, written for i = 1,2, 3, with a specific rule for gauge 
transformations. Similarly, e is locally a 1-form over M transforming through the adjoint representation. 
If g is map from M to SU{2), the corresponding gauge transformation reads A i— )• Ad{g)A + gdg~^ and 
B I—)- K(l{g)B. The adjoint action is in matrix form: Ad(g)X = gXg~^ . To write the action, note that su(2) 
is equipped with a non-degenerate bilinear invariant form, which enables to raise and lower indices via the 
Kronecker delta 5ij. The curvature of A is the su(2)-valued 2-form F{A) = dA + and the action is 



•5(e,^) = 7^ / e'^F{Ay 5ij . (3) 

IM 
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In the following, we will set AttG = 1. 

The equations of motion are the torsion free condition, e = de + [^4, e] = 0, and the flatness condition, 
-F = 0. It can be shown that all solutions are locally equivalent to the trivial one, A = Q and e = 0. This 
fact is due to an additional symmetry of the action, 

A^ A , e t-^ e + dAV 1 (4) 

for any function r/. This symmetry makes possible to locally gauge the cotriad to zero. 

The Hamiltonian analysis reproduces the covariant picture, while preparing the road to quantization a 
la Dirac. To see how this happens, assume M is of the form S x R, where S is a smooth two-dimensional 
manifold of arbitrary topology. Choose arbitrary coordinates = (x^,x^) over the canonical surface S, 
together with a 'time' coordinate t on IR. Following this splitting, the action admits the decomposition 

S{e,A) = j dtj d^x [e'^H.j {^adt^,) + 6ijAi [e'^' D,el) + 6,,el [e^' F^,)] , (5) 
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with e"'^ = e"''*. The operator D is the covariant derivative for the restriction of vl to S. The interpretation is 
straightforward: the canonical pairs are formed by the components of the connection and their momenta 



Ef = Sije""^ ei, the fundamental bracket being 



{A^,{x),E'^{y)]=5lS^^5^^\x-y) . (6) 

Furthermore, the time components et and At appear as Lagrange multipliers, enforcing the Hamiltonian to 
be a combination of constraints 

DaEt = 0, F^, = 0. (7) 

The first one is the GauB law, which generates the SU(2) gauge transformations on the canonical variables. 
The second equation, F^^ = 0, constrains the connection to be flat and generates the spatial part of the 
special symmetry (jj]) via the brackets. Together, these constraints form an ISU(2) algebra. 

When the frame field is non-degenerate, it is possible to extract from the multiplier et the usual lapse 
and shift variables, so as to split the constraints F^^ = into a part generating diffeomorphisms of S and 
another part called the scalar, or (somehow loosely) Hamiltonian constraint, (s^. To this aim, we introduce 
the normal density vector 



n 



'-e'^^eabE<^El = {e^xe2)\ (8) 



using the vector notation for internal indices. Also define the curvature density vector, F'^ = e"'^F^f^. It 
turns out that the norm of n is exactly the determinant of the 2- metric on T, induced by the co-triad: 
|np = det(^g'). Then, since n*£'f = 0, we can write 

et.F = N''Va + — , ^ NH , (9) 

where A'^'^ and N are respectively the shift and the lapse, defined as 

= ^ \ , e\fc Et ei N= ^ ^ et ■ n . (10) 

det(25) * * TditM 

This way the vector constraint Va, which imposes the diffeomorphism invariance over S, and the scalar 
constraint if, which is simply the projection of F onto the normal n, appear 

Va = E\F:, = ea-F , (11) 
H = etFl,E'^El = ^^.F . (12) 

At this stage, an important remark is in order. In contrast with the proposal of 0], we will consider a 
Hamiltonian H of density weight 2 rather than one. The reason for this choice is that we will perform the 
most naive regularization and quantization, in which E'^ and F'^^^ are easily treated, but not non-analytic 
factors like —1===. 
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B. Loop quantization 

We now proceed to the loop quantization on the unreduced phase space, and we wih then study the 
reduction by the constraints at the quantum level. First, the construction of the space of states in the 
connection polarization is based on cylindrical functions, which probe the connection only through a finite 
number of variables: the holonomies Ue{A) of the connection A along some set of edges (e). Given a graph 
r, with E edges and V vertices, and a function / over (SU(2))-^, we form the cylindrical function iprj as 

^rj{A) = f{UeM),---,Ue,{A)) . (13) 

We are moreover interested in SU(2) gauge invariant states. Gauge transformations act on holonomies only 
on their endpoints. If /i is a map from S to SU(2), then: Ue{A^) = h{t{e)) Ue{A) h{s{e))~^, with t{e),s{e) 
being respectively the source and target points of the path e. When focusing on a single graph T, this 
reduces gauge transformations to an action of SU(2)^ on the set of cylindrical functions over F. So from 
any function / over SU(2)^, one gets an invariant function by averaging over the SU(2)^ action. 

The algebra of such functions has a natural inner product which comes from the Haar measure on 
{S\J (2))^, dfir = Hedge, giving 

{(l>r,h\iprj) = Y[dge h{gi, ...^qe) figi, ■■■,9e) ■ (14) 

e=l 

Completing this space for the corresponding norm leads to the Hilbert space L^(SU(2)^/ SU(2)^, d/xr). 

To get the full kinematical Hilbert space of loop quantum gravity, one has to consider a sort of union 
of cylindrical functions for all possible graphs. This results in a Hilbert space which can be understood as 
L^(^, dfiAh), the space of square-integrable functional over A. An element of ^ is a generalized connection, 
which assigns a holonomy to any path, while the Ashtekar-Lewandowski measure ^al, arising from the set 
of measures fir , makes cylindrical functionals on different graphs orthogonal (see for instance 0, Q ) . 

The spin network basis is an orthogonal basis of -L^(SU(2)^/ SU(2)^, d//r) • The first idea is to expand 
any element of L^(SU(2)^, d^r) onto the matrix elements of the irreducible representations of SU(2)^, 
labelled by E spins je G N/2. Then, rotation invariance at vertices imposes to contract the matrix elements 
of edges meeting at a vertex with intertwiners l^. If w is a node of F with ingoing links (ein), and outgoing 
links (cout), then an intertwiner is a map: ^ei^T-ij^ — >■ •S'eout^ie) (where Hj is the carrier space of the spin j 
representation) which commutes with the group action. Thus, to form the spin network basis on we need 
to assign spins to edges and a basis of intertwiners at each node. For a trivalent vertex with ingoing edges for 
example, there is a single intertwiner /-jijajs • T~ij^®l-Lj2®T~Lj.^ i— t- C up to normalization, its components in the 
standard magnetic number basis being the Wigner 3jm-symbol {ji^nii; j2,m2; j^,m2,\ij^j2j-j,\Q) = ma)' 
Hence, spin network states are the following SU(2) ^-invariant functions 

E V 

Sr^'^'Hai, ■ ■ ■ ,9e) = ^ JJ(je,?7ie|fi(e|je,ne) JJ((g)em jeimekt;! (Xieout je,ne) , (15) 
{me,ne}e=l ""=1 

which form an orthogonal set. This enables a transform between gauge invariant functions on F and functions 
over the colorings {je,iv)- 

Let us now look at the algebra of operators. Quite clearly, spin network functions form an algebra, so 
that a gauge-invariant holonomy operator can be defined as acting on T-Lr by simple multiplication. In order 
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FIG. 1. 



to form quantities depending on the triad Ef, it is natural to smear it along curves, thus defining the fluxes 

X,j = jEt r eabdx\ (16) 

for a path c and a function /. This is the natural choice given the chosen smearing of the connection so 
that the Poisson bracket closes on the space of cylindrical functions (see typically [13]). It can be shown 
that the bracket {X(.j,il^T,h{^} receives contributions only from the crossings of c with F, and that there 
the flux act as a left or right invariant vector field on the function h. 

This framework allows to define different kinematical observables such as the lengths of curves, and areas, 
like in [2]. As we will need some of these, we will proceed to a simplified derivation based on the phase space 
of loop quantum gravity restricted to a single graph, section [TVl 



C. The flat state on the boundary of a tetrahedron 

For concreteness of the discussion, we consider the 2-sphere as our canonical surface, triangulated with 
four triangles which are glued like on the boundary of a tetrahedron. Our computation actually works 
without amendments for any 3-valent vertex of a triangulation. It can also be extended to vertices of higher 
valence (by duality, to cycles of spin network graphs containing an arbitrary number of links j^]- 

The kinematical Hilbert space of gauge invariant states is spanned by spin network states on the oriented 
graph dual to the triangulation, which is also a tetrahedral graph. The states can be seen as functions 
over the six SU(2) elements on the dual edges and invariant under translation acting at the four dual vertices. 
Let us write explicitly the spin network function, labelled by six spins (ji, . . . , jg), following the figure [H 



{gi,...,ge)= ^ 



ai,...,as 
bi,...,be 



31 

h 



32 
-02 



J3 

-as 



31 35 36 
-ai 65 -qq 



33 3a 35 \ (32 36 3a 



■ (17) 



^ Though, the recurrence relations one gets in these generalized situations are more intricate, and require further discussion 
beyond the scope of the present note (see 
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The range of summation is —je < Oe,6e < je for each hnk. In the sum, each Hnk e carries two magnetic 
indices, be on the source vertex which is contracted with the right of the holonomy g^, and on the target 
vertex and contracted with the left of ge- 

Given the trivial topology we consider, there is only one physical state, which is not normalizable. It 
formally corresponds to the product of Dirac deltas for the holonomies around the independent cycles of 
the graph. Cycles can be defined as the boundaries of the faces dual to the vertices of the triangulation. 
The product has to be considered over independent cycles, which in the present case is easy to understand: 
asking for the holonomies around three faces to be trivial enforces the same condition on the remaining, 
fourth face. So we choose any three cycles like (see pictured]) 

V'physC^i, • • • ,96) = ^{§4^ 9596)^(92 '^91^96)^92^ 9-m) ■ (18) 

To get its expansion on the spin network basis, we take the inner product with a spin network function 
sli;\gi,...,g6), 

V'phys(ji, • • • , je) = [Hdge 4et^(9i, ■■■,96) V'phys(5i, ... ,96) = 1^.^ f] ■ (19) 

e=l 

The quantity into brackets is known as the Wigner 6j-symbol. This result simply comes from the fact that the 
physical state enforces the evaluation of the integral on the flat connections (up to gauge transformations), 
while there: ^{^^(l) = {6je}- 

The 6j-symbol is known to satisfy the Biedenharn-Elliott, or pentagon, identity 



where S is the sum of the nine spins. It is key to the Ponzano-Regge state-sum model, since it corresponds 
to stating its invariance under a change of the spacetime triangulation by a 2-3 Pachner mov^. 

Specializing (|20p to the case / = 1, one is free to choose g = d and h = c, while the sum on the right 
hand side reduces to only three terms, for / = A: — 1, fc. A; -|- 1. This leads to the second order recurrence 
relation 

A,^(n) 1^*; ' f f I + MM If f f] + A_,(M - ' f f I = . (21) 

The coefficients are 

Ao(ji) = (-iy^+^'^+^^l^^ Ij + (_i)2ii+j2+j3+j5+i6(2ji + l) l^} -^i U U , (22) 

IJ6 J6 J4j [j2 J2 JZ] Ue J6 J5j 

A±i{h) = (-I)2ji+i2+i3+i5+i6+i(2(ji ± 1) + 1) {-^i ^} U j-^'i ^} U . (23) 



They can obviously be evaluated explicitly (but the above expressions will be more useful to us), and the 
recurrence relation can be solved, starting from an initial condition to determine all values of the 6j-symbol. 



^ To prove topological invariance, one also needs to check the invariance under the 1-4 Pachner move. It can be done using 
again the Biedenharn-Elliott identity and some orthogonality properties of the 6j-symbol. However, the latter invariance is 
only formal since it results in a divergent sum, [sl. ITst. 
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An interesting feature is that though it is a second order difference equation (involving three consecutive 
neighbors), the solution is fully determined by a single initial condition. Looking at the lowest possible value 
of the spin ji, that is: = max(|j2 — js], |j5 — jel), it turns out that the lowering coefficient 
evaluated on this spin is zero. Thus, the recurrence can be implemented starting from the initial value on 
this lowest spir|l. 

So, the recurrence relation and the Biedenharn-Elliott identity 

• encode the invariance properties of the Ponzano-Regge model, or equivalently they are responsible for 
its special symmetries, 

• and fully determine the 6j-symbol which is the physical state in the spin network basis for the trian- 
gulated 2-sphere. 

To make these two aspects talk to each other, we would like to derive the recurrence relation as an equation 
encoding the symmetries in the context of loop quantum gravity, in the form of ([2]) . 

III. THE PROPOSAL FOR A NEW HAMILTONIAN 
A. The discrete geometries of loop quantum gravity on a single graph 

When looking at a graph, the loop quantization, which is based on specific smearings of the fields, can 
be reproduced starting from a classical phase space associated to the graph. The geometric interpretation 
of this phase space will be of particular importance for us. 

Let us restrict attention to a single oriented graph, typically the dual to a triangulation of S. The 
phase space inherited from the smearings which are used in loop quantum gravity is the cotangent bundle 
T*(SU(2))^. More precisely, the pairs of canonical variables (^^,£'f) are replaced with a pair {ge,Xe) G 
SU(2) X su(2) on each edge of F. The group element is the parallel transport operator along e from its 
source vertex to its target vertex, while is the flux variable for an edge dual to e. We chose to act as 
a right derivative, which means that we have the brackets 

{Xl,Xi} = i\Xl , {X:,5e} =9eT\ (25) 

where the matrices (Tj)i=i^2,3 are anti-hermitian matrice^ satisfying [t^t-^] = ^t^ . All other brackets 
vanish. The right derivatives are obtained by transporting Xf. to the target vertex of e, 

Xe = Ad(5e)^e, together with {5C^,Xi\ = -e'\X'l, {Xl,g,} = . (26) 

Here Ad is the standard adjoint action of the group on its algebra, and in particular, it is here an action of 
holonomies on fluxes. 

^ This property is important with respect to the asymptotic hmit. There, it is known [lol. [igl. [20I] that the 6j-symbol oscillates 
like the cosine of the Regge action SaQ'e) = ^ZeO^ + for a tetrahedron with lengths (je + |), 

fjl h d-i\ ^ 1 1 /p»(SR(je) + f) I p-»(SR(je) + f )\ (■24) 

Having a second order equation, it may have been expected to be solved in the asymptotics by both the positive and negative 
exponentials of the action. But the model is fully independent from the orientation: there is only one required initial condition, 
because of the vanishing of leading in the large spin regime to the cosine of the action. 

They read: = — f cr*, in terms of the Pauli matrices (ct*). 
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In this framework, the Gaufi law becomes the sum of the fluxes meeting at a vertex u of F 



(27) 



if all edges are outgoing. Using this equation, we can get a nice interpretation of the flux variables as describing 
the discrete intrinsic geometry of the graph dual to T. For instance if F is dual to a triangulation of S, it 
is natural to see the variables as 3- vectors describing the edges of the triangulation in R^, so that the 
above equation is nothing but the closure of each triangle. 

One can then extract the intrinsic geometry as follows. The length of an edge e* of the triangulation 
is the norm of Xg, for e being the dual link to e*, ie = \Xe\. The dihedral angle (pf,^' between two edges 
(e*, e *) of a triangle is contained in the dot product of Xe,Xf,i, where the dual edges e, e' meet at a vertex. 



where e^e' is 1 if e and e' are both outgoing or ingoing, and —1 else. 

It is also possible to introduce some notions of 3d dihedral angles, i.e. between two triangles sharing an 
edge of the triangulation. Since it will be a central point in our discussion on the dynamics, we postpone 
its presentation to the next section. 



The usual flatness constraint on a triangulation, inspired from Regge calculus, is to consider curvature 
concentrated around {d — 2)-simplices, i.e. here vertices of the triangulation. In the dual picture this 
corresponds to looking at the parallel transport operators around the faces of F, and ask them to be trivial, 
so that 



Together with the Gaufi law (j27p . they form a system of first class constraints, in which the above flatness 
constraint generates translations in of the vertices of the triangulation. This is the usual starting point 
when one wants to reduce the phase space before quantization, since it is easy to see that the theory has 
actually only a finite number of degrees of freedom (given in terms of the topology of S). 

This is also the form of the constraint which has been mainly used to reduce after loop quantization and 
make the link with spin foams, like in [^,0]. From the continuum LQG point of view, the underlying idea 
is to define the projector onto flat connections through its matrix elements between any two spin network 
states. This is done using a lattice regularization which is adapted to the spin network graphs one wants 
to consider, like a 3-valent graph F, dual to a triangulation and which contains the spin network graphs. 
Then, one has to define the flatness constraint on plaquettes of F, ()29p for instance, and this results in 
the projection of any spin network whose graph is contained in F onto its restriction to flat connections. 
When evaluating physical inner products, the topological nature of the moduli space of flat connections 
guarantees the independence with respect to the triangulation used to regularize the constraints. Thus, it 
can be checked that the regulator can be removed [^. We will use implicitly this reasoning and hence simply 
work on a fixed graph F, dual to a triangulation. Since we decide not to use (f29]) . we have to see if our 



cos 4>ee' 




(28) 



B. Construction of the new Hamiltonian 



(29) 
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approach makes us able to recover the projector onto physical states on T. This will be done partially in 
the section [TVl for trivalent faces. 

In contrast to the form (j29p . we are here mostly interested in a form of the constraints which could 
be closer to what we need for 4d gravity, like in Q], and at the same time make contact with recurrence 
relations which are a way to encode the symmetries at the quantum level in spin foams [11]. Our proposal 
is to quantize the following operator in the spin network basis 

H,j = Xe,-X,,-X,,-Ad{gf){Xe,)= Yl K{^iJ-Ri9fh)xi. (30) 

i,j=x,y,z 

and show that the corresponding constraint equation is exactly the recursion relation on the 6j-symbol when 
the triangulation is the boundary of a tetrahedron. Here 61,62 meet at the vertex v of F, and are in the 
boundary of the face /. In the first line, we have used the adjoint action of the group on its algebra, and in 
the second line, we have emphasized that it is equivalent to the vector (spin 1) representation, given in the 
Cartesian basis by the usual rotation matrices. 

First we give a nice geometric interpretation of this quantity in terms of discrete geometries, which in our 



view strengthens and makes more precise some of the ideas of [2l|, |22] which were proposed at the classical 
level to describe a flat 4-simplex from its boundary. 

There are two basic ideas which support our choice for H^j. The first key idea is to represent the scalar 
constraint H, in the following way. As can be seen in (jl2p . the constraint is just the contraction of the 
curvature with the normal n, ([8]), to the canonical surface. Given the choice of smearings in LQG, the 
curvature on a spin network graph is naturally regularized as the product of the holonomies along the cycles 
of the graph, where a cycle is defined as the boundary of a face delimited by T on the canonical surface. In 
practice the curvature tensor is replaced with 

e\F^, 5^^-{Rigf)f. (31) 

Equivalently, the curvature is the effect of parallel transport around a vertex of the dual graph to T. But 
then, it is quite unclear what a normal to this vertex means. So we choose a face touching this vertex and 
consider the normal to this face. Back to the description on the spin network graph, the latter is dual to 
a vertex of F, and choosing a normal Nt is equivalent to choosing two flux variables X^^ and Xe2 meeting 
there and both lying in the boundary of a common face, with: Nt = X^^ x . So the natural choice is to 
contract the regularized curvature along a cycle of F with X^,^ and X^^-, where 61, 62 are edges of the given 
cycle c = df meeting at a vertex v, 

EtE] K^Xi^. (32) 

Also, we choose v as the base point to define gf and take its orientation so that it transports first along 62 
and ends with 61. Then 

EtEU'\F^^, H^j. (33) 



One may ask then about the splitting of the constraint (j29p into scalar and diffeomorphism constraints 
and the precise status of our proposal in this respect. Remember this splitting comes from projecting the 
curvature onto the normal n to S on the one hand, (I12p . and onto the directions tangent to S on the other 
hand, (jlip . Here assuming we are dealing with a triangular face /, we have three different H^j, one for 
each vertex. Each of them corresponds to choosing the normal to one of the three triangles meeting at the 
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node dual to /. So for generic X^, the whole set of H^j for a given face takes into account the projection 
of the holonomy around / in the three independent directions of R,'^. For this reason, there is no need to 
introduce diffeomorphism constraints in this framework. 

A crucial fact to make our assertion complete is to check that the algebra of the set of {Hyj) on F closes. 
Here comes the second key idea of our construction. If one wants to impose the flatness constraint as a 
Wheeler-DeWitt equation, one needs to find three (real) independent constraints on the matrix elements of 
the holonomy gf. To do it in a gauge invariant way, one has to project the matrix R{gf) onto the fluxes, 
and this is how H^ j is constructed. 

Consider a cycle with three edges, say ei, 62, eg like in the figure[Tl with holonomy gj = gQ^gig2- The three 
functions H^j on this cycle are actually the three constraints obtained by contracting (id —R{gf)) with pairs 
formed among three vectors defined at the node v where 62, eg meet, which X2,Xq and Xi{v), where Xi{v) = 
Ad{g2^)Xi. We denote the Hamiltonians by H26, H21, Hiq. For example, H21 = X2 ■ {id—R{gf))Xi{v). 

In this sense, this is just looking at the matrix elements of R{gf ) in a non-orthogonal basis. For a generic 
triangulation, there is one constraint for each node of a cycle. Thus there are the appropriate number 
of constraints to enforce R{gf) = id G S0(3), when the flux variables around the cycle span the three 
dimensions of the algebra. 

To complete the argument, we first show that the three H^j around / are good coordinates in the 
neighborhood of the solutioi|l/?(5j) = id. Writing g{p) = — p? \ + ip ■ a € SU(2), we have 

X ■ X' - X ■ R{g)X' = -2V1 {X X X') ■ p + 2[p^ {X ■ X') - {X-p){X' -p)] . (34) 

From this one finds \dH2Q A dH2i A d^^i6|_R(g/)=id = 8| det{X2, Xq, Xi{v))\'^ , where the differential is taken 
with respect to pf. Hence the system is generically of rank three, and this makes clear why we need the 
three vectors to span a basis of . 

Finally, all Poisson brackets weakly vanish. Consider the bracket between H12 and i^ig which are in the 
same face. It reads 

{Hu,Hie} = {Xq X X2) ■ [id-R{gf)]X^{v) + {X2 + X^) ■ {Xi{v) x R{gf)Xi{v)) , (35) 

which is zero when R{gf) = id. Consequently it can be re-expressed as a linear combination of H2Q, H2i,Hiq. 

On the boundary of the tetrahedron (figure[TD, the other non-trivial brackets are {H12, H23}, {H12, H^^}, 
{Hi2,H24} and {H12, Hq^} (and the other are deduced by symmetry, or identically zero). It is straightfor- 
ward to check that they can be written with the quantities {id —R{gf)) and Xg x R{gf)Xe. 

C. Its geometric meaning 

Physically, what we expect from the classical flatness constraint F^^ = is that S with its geometry 
can be locally embedded into flat 3-space. From the point of view of the geometry of a triangulation, it 
means that the (3d) dihedral angles Qtit2 between two adjacent triangles ti,t2 should be given as a standard 
function of the dihedral 2d angles (0ee')- ^ typical situation is depicted at figure [21 for a trivalent node of 
the triangulation. The triangles ti,t2 are dual to the vertices vi,V2 of the spin network graph F, which are 



^ Notice that the solutions of the equation X2 ■ — X2 ■ Ad{g)Xe, = are g = exp{t2X2) exp{rj{X2 x Xa)) exp{tfiX(i) , where 
ti,t2 are arbitrary. But rj admits only a finite number of values, since there is a finite number of S0(3) rotations with axis 
X2 X which solve the constraint. It seems that non-zero ri however only survives if X2 x Xe is proportional to the third 
vector Ad{g2^)Xi. 
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linked by the edge e. The classical angle between the two triangles ti,i2 can be computed from the three 
2d angles around the node c* , which is dual to the cycle c, using the formula 

r> V\ COS(/>eie2 - COS(/>eie COS(/>e2e .„„^ 

COS etit2 [X, X) = — — . (36) 

We have written it as a function of the flux variables Xe,Xe since all angles (/>ee' can be evaluated from 
them, without writing down holonomies. In addition, when the closure relation (j27p holds classically, the 
angles (pee' are determined as functions of the lengths {£e)i and hence so are the angles Qtit2- 

Parallelly, using the interpretation of the flux variables and the holonomies, there is a natural notion of 
dihedral angles between triangles, defined as follows. Since the fluxes represent the directions of the edges 
of the triangulation, the normal Nt to a triangle is given by the wedge product of two of them. Following 
picture [21 A''^*^ = e^-^^Xi X^^, and Nl^ = e^-^^Xi^ X^. It suggests to look at the dot product of these normals 
as being the cosine of the dihedral angle between ti and t2- But for this, it is necessary to transport them 
in a common frame, say Nt^ along the edge e to the vertex vi. So we define the angle 0t^t2 as 

n (Y Y ^ti ■ Ad(ge)jVf2 

COS 9t,t2 [X, X, g) = I I — . (37) 

This is a function of the holonomy since one has to compare the local embedding of ti to that of t2. The 
result can be partially evaluated in terms of the 2d dihedral angles, because it is a dot product of two vector 
products in I^? 

a fv V \ ^l-^'^{9e)X2 - COS(/>eie COS(/)e2e .ooN 

cos et^t2 {X, X, g) = — — . (38) 

Sm^eie Sm^eae 

Quite clearly, the dependance on the holonomies cannot be dropped off. Obviously, one would like to 
compare this notion of dihedral angle to the standard formula (I36p for the angle Qtit2- It follows from the 
definition X^.^ = Ad(5'ei)^i (and the same for e^) that the difference between 9t^t2 and 6*1*2 is precisely 



H^j = sin4>e^e sin0e2e (cosGtita - cos 6**1*2) ■ (39) 
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We will not go further about the classical aspects, though a full classical analysis would certainly be 
worth. It could be done either with the variables {g^X) or with the gauge invariant variables {1,6). In this 
respect, it is interesting to note (see 0, [i^]) that the angle Of. = 9t-^t2 {X,X,g) associated to the edge e is 
in fact the canonical variable dual to the edge length, 

{/e,ee'} =5e,e' • (40) 



IV. A WHEELER-DEWITT EQUATION FOR THE FLAT STATE ON THE BOUNDARY OF A 

TETRAHEDRON 

Now we look at the quantization of H^j on the boundary of a tetrahedron, see figure [H for v the vertex 
where j2, J4,i6 meet and c = df the cycle {jij2j&), 

Hv=(e2nea)J=ijij2j6) = X2-Xq - X2 ■ Ad{g2^g^^ge){XQ) . (41) 

This is done is two parts: first the quantization of the operator {X2 ■ Xq), which is basically the cosine of the 
angle between the edges 62,66, and then the same quantity computed after parallel transport of Xq around 
the face. In both cases, we use the space of gauge-invariant states Tir, spanned by spin network functions. 



A. Quantization of the cosine of the 2d dihedral angle 

The observables of the intrinsic geometry, lengths and 2d angles, are easily defined on spin networks. 
Firslll, 

{Xl 4'"'"^) (51, ...,gE) = -{Xl, {Xei,4'"^''^(5i, ...,gE)}}= jeiJe + l) 4'^''"^(5i> ■ ■ ■ , 9e) , (42) 

for the lengths, diagonalized in this basis with a discrete spectrum. 

As for the dot product of different Xe meeting at a node, it can evaluated as follows. We consider the 
edges 62,66, meeting at a 3-valent node with 64, and assume without loss of generality that they are all 
outgoing. Then, the product X^X^i acts by inserting r^^ iS) Te^i in the spin network, thus modifying the 
intertwiner. This corresponds to adding a link in the spin 1 representation between the two edges, on the 
right of the group elements 562)566- In equations, we have 



{X2-Xg) V (-^^ "'^ ■^^] 5e2|j2,a2) Leslie, ae) '^c/e4li4, 04) 

a2,ai,a(j 

(???) 9e2T'\j2,a2) geenlh^ae) (g) ge4\ji,a4) 
^ \a2 ttQ a^J 



02,04,05 



To extract from this formula some Wigner recoupling coefficients (a 6j-symbol actually), we need a precise 
calculation of the so-called grasping r'^ Te^i- 

It is convenient to go to a different basis of the Lie algebra, in which the flux variables X^ have magnetic 
indices m = —1, 0, 1 instead of the Cartesian indices {i = x, y, z), and with generators Lm- A way to choose 



The scalar product is taken to be X ■ y = —X^Yi because we are using anti-hermitian generators. 
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well adapted generators is the following. We have used so far the adjoint action of holonomies on fluxes, 
which is equivalent to the spin 1 representation. But at the quantum level, the matrix elements of the adjoint 
action act by multiplication on cylindrical functions. Because we built the spin network functions in terms 
of Wigner matrices, it is natural to choose a basis where the matrix elements of the adjoint representation 
are the Wigner matrices with spin 1. 

We can test this request on the parallel transport relation (j26p at the quantum level. In the Cartesian 



basis, we have: (Ad((7e)Xe)* = R{gey jXi, where R{g) is the standard rotation matrix. When acting on 
matrix elements (of the spin network basis), that becomes 

RideTj Xi (j, 6| ge |j, a)=Y, RideTj {j, b\ ge T'j I j, a) = {j, b\ ge \j, a) = x\ {j, b\ ge | j, a) . (44) 



is 



We have used the well-known intertwining property: R[gY ^t^ = g T^g. Thus, the classical relation 
correctly implemented on wave functions. We now ask for the same relation to hold in the new basis (^^772) 
and with the Wigner matrix of spin 1 instead of the rotation matrix R{g), 

{j,b\ge^ Lrrige\j,a) = ^ {1, m\ g \1, n) {j , b\ \j , o) . (45) 

n=-l,0,l 

That shows that the matrix-valued vector L in the spin j is actually an intertwiner T-Lj T-Lj* Hi — ?■ C, 
where Tij* is the dual representation. Hence, the element {j, b\ \ j, a) must be proportional to the Wigner 

3mj-symbol (— 1)-'"'^ ( ^ -air choose the normalizatior0 like 



{j, b\ Lm \j, a) = (-1)2^+1 N, i-iy-'^ , (47) 

the prefactor being 

Nj = ^j{i + l)dj , with dj=2j + l. (48) 

The insertion of the generators in (j43p can be done using 

+1 

Ti®T'= {-If-"^ Lm®L_^. (49) 

m=— 1 

By plugging the above formulae into (|43|) . one gets the following modified intertwiner in terms of 3jm-symbols 
(j2,a2; j6,a6;j4,a4|2:j2j6i4|0) = i!"^) 02,021 -^^m |j2,a2) (i6,a6| L_m (-1)^"™ |j6,a6> 

^ ' \ai a^J ^ ' \m -02 a^j \m -a^ a^J ^ ' 



Notice that these generators are quite unusual since 



Lo\j,a)=a\j,a}, and L± | j, a) = T\/ ^^^^^^^^r^^^^ b> T 1) • (46) 
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Prom the invariance under rotations on the indices (02,05,04), it fohows that this intertwiner must be pro- 

' 32 je 34 
aL clL a4 



portional to the standard 3jm-symbol { ''J ''J ) ■ '^^^ proportionahty coefficient is obtained by contracting 



^hieh "with this 3jm-symbol, and is found to be a Wigner 6j-symboil, so that 

((iTle) sli^){g,,...,g,) = {-ir+^-+^-N,,N,, \f f U 4'^^(5i, • • • ,56) • (52) 

The fact that we find a spin 1 in the symbol is simply due to the fact that the indices contracted in the 
grasping ()49p live on the vector, i.e. spin 1 representation. 

The above expression will be very useful for us. Still, it can be simplified thanks to the known value of the 
above 6j-symbol, and observing that among the prefactors we find the eigenvalues of the length operators 
for 62, cq. This naturally leads to 

/ —T- / N [hih + 1) + Mh + 1) - jiUi + 1)] {j,}, . 

(cos 06266 slit )i9i,---,9(i) = n /■ r- -lA • / ■ 1N 4 (5i,---,56) ■ (53) 

2VJ2(J2 + 1) Jeije + 1) 

This exactly reproduces the classical expression of the dihedral angles of a triangle, with quantized lengths 
ieije) = \/ je{je + 1)- Notice that a similar expression, with a prefactor going to 1 in the large length limit, 
was found though in a completely different setting based on discretizing the path integral on a triangulation 
of spacetime [25| . 

At this stage, we are able to notice that the part {X2 • Xq) of H^j produces a 6j-symbol and a phase 
factor which are exactly those entering the first term of the coefficient ylo(ji), (|22p . which appears in the 
recurrence relation on the 6j-symbol. 

B. Quantization of the cosine of the 2d dihedral angle after parallel transport 

The next step is to compute the action of the same dot product after the insertion of the holonomy along 
the cycle. We first rewrite this quantity using the left derivatives X^, 

X2 ■ Ad{g^'g^'ge){Xe) = X2 ■ Ad{g^'){Xe) = ^ 1^ {6i, - R{g^%) xi . (54) 

Then, at the quantum level, X2 g will insert a generator r* on the left of the holonomies g'2,6, while R{g^^) 

will act as a multiplication. Let us focus on the part of the spin network function s^g^ involving the matrix 
elements (j'l, oi|(7i|ji, 61). Then, the operator X2 ■ Ad{g^^){XQ) acting on 



^ \-ai 05 -oe/ \oi -02 -03/ 

"■iM ^ ' ^ ' (55) 

«2,3,5,6 

b5,a5) ^ (isjOsI O (J6,a6b6 ® {j2,a2\g2 , 



The result comes via the following formula 

jl h h\ ^y^(^_^^^Ul{j^-a,) fjl h J3\ / js J4 j5\ f jl js j6\ f J2 jd ji\ 

jl is ie J ^ \ai a-2 as J \-a3 —04 as / \— tsi -as ffle / \-a2 — ae 04 J ' 

which is a standard definition of the 6j-symbol. 
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gives the following expression 

f-^ \—ai as —gqI \oi —02 —03/ 

ai,oi,m,n 
12,3,5,6 

Ii5,a5> <^ Os^OsI (g) (j6,a6|-^^-n(-l)^~"56 {j2,a2\Lm92 ■ 

(56) 

Here we have re-expressed the insertions of generators in the basis {Lm), which transforms with the Wigner 
matrix of spin 1, and we have used the fact that (— 1)™'~"(1, — m|(7^^|l, — n) = {l,n\gi\l,m) . 
The product of two matrix elements of gi is expanded onto irreducible representations via 

ji+i 

(l,n|5(i|l,m)(ji,ai|5i|ji,6i) = (-l)^^i ^ dj^{Ji,Ai\gi\Ji,Bi) 

Ji=ji-1 Ai,Bi 

1 jl Jl \ I i\Ji-Ai ( 1 Jl Jl \ ( ^\l-m+ji-hi 



n ai -AiJ ^ ^ \-m -bi Bi) ^ ' ^ \ 1 

where Ai = n + ai, Bi = m + 61, while the action of the generators produces 

(i6,a6|(-l)^-"L_„ = (-l)2.«+iiV^.^^(-l)i-+^a-4 ( 1 ^6 ^^.^^^,| _ (59) 

The contraction of the magnetic indices between the above formulae results in the addition of a link in the 
representation of spin 1 between the links carrying Ji and j2, and similarly between Ji and jg. It is then easy, 
using the same method as that used in the calculation of the dihedral angle, to extract 6j-symbols at the 
two vertices of ei, in order to re-express the function in the spin network basis. After these manipulations, 
one gets 

[X2-M{g^^g,){X^)) 

Ji+l r r ■ r . ■ (60) 



Jl=jl-1 



h 


Jl 


is} 


y 


Jl 


Js} 


{J2 


J2 




[J6 


J6 





'tet 



The prefactors, Nj^Njg, are the same as those arising in the action of X2 ■ Xq. In addition, one can recognize 
the coefficients ^±i(ji), (j23|) . in the terms Ji = ji±l, and the second part of the coefficient Ao(ji), ([22]) . for 
the term Ji = ji. The technical point which requires to be careful is the sign of these coefficients. Finally, 
one obtains 



where the coefficients Am{j) are given by the equation ([22 



^OlJejStet +^-l(jejStet + ^+UJeJ ^tet 



(61) 
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The last step is to derive the constraint equation ([2]) for a hnear combination of spin network states 

6 

V'(5l,...,56) = [n(2ie + l)]V'(jl,...,i6)4et^(5l,---,56) ■ (62) 

A brief examination of the coefficients ^±i(ji) reveals that: (iji=pi^±i(ji =F 1) = '^ii^Ti(ii)- This way, the 
equation Hvj'4){gi, . . . ^qq) = exactly becomes the following difference equation of the second order 

A+i{je) ipih + 1, . . . , je) + Aoije) V(ii, . . . Jg) + ^-i(je) " 1, • • • , je) = . (63) 

This is obviously the recurrence relation which defines the 6j-symbol from the Biedenharn-Elliott identity. 
We emphasize here that we have produced the exact equation, and not only a large spin approximation. 

C. Equivalence with flatness on trivalent plaquettes 

We now want to go further that a single tetrahedron, by showing how the Hamiltonian we propose enables 
to recover the projector onto the trivial holonomy sector for trivalent faces. This process can be thought 
of as a Pachner move 3-1 which reduces the triangulation by removing a dual face. Let us explain how 
our Hamiltonian implements it. In the dual picture to the triangulation, the initial configuration is a cycle 
bounding a face with three links. Thus we have three constraints H^j, that are good coordinates in the 
neighbourhood of R{gf) = 1. At the quantum level we get three difference equations to impose the three 
constraints on the holonomy gj around the cycle. 

We take the same notation as before, denoting ei, 62, eg the links along the cycle, and 63, 64, 65 the links 
which connect to the cycle. This is like in Figured! except that the links 63, 64, 65 do not necessarily meet. A 
key feature of the difference equation, which can be observed in the tetrahedral case, is that it only depends 
on the spins on the cycle and on those which connect to the cycle. It is independent of the fact that the 
latter (63,64,65 here) meet. 

Using the three difference equations, we first want to show that the physical state factorizes on the cycle. 
In the spin network basis, the state is a function of the spins. In this configuration, its dependence on the 
spins ji,j2,j6 actually completely factorizes, 

^phys(je) = (-1)'^'^ If f fl ^ife), (64) 

where the spins jg are the spins of all links of the graph except 61, 62, 66- This is known in the 3d gravity case 
from the projector onto physical states [6], and work is in progress towards the precise loop quantization 
of higher dimensional BF theory [26]. Such a factorization property is actually very natural. Indeed, a flat 
connection induces on the boundary of a face with three links holonomies which are necessarily trivial up to 
gauge. Then, it is known that spin networks evaluated on the identity satisfy such factorizations, and the 



case of cycles with more links is also known [111] • 



The idea to extract the factorization is that the recursion relations hold on ii,j2)j6- We can thus 
implement them from ip{ji, j2, je, j'e) so as to reach an initial state with ji = 0,^2 = js and je = js. The 
factorization (I64p is then proved with the initial conditior 1^ V'(ie)|ji=o = (—1)^-'^ | | 4'ije)- 



^ It will be convenient in the following to explicitly factor a phase (—1)'^-'^. The basic reason is that in our conventions 63 is 

oriented inwards the cycle while 64,65 are outwards. 
1° Notice that { ^° « « } = 



IS never zero. 
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Now we are ready to present the consequences of the factorization ()64p in the group representation of 
wave functions, in which we also get a factorization of the cycle, 

V'phys(5'l,5'2,5'3,5'4,fi'5,S'6, • • • ) = ^{92^ Oi^ Og) 4'{g2^ 93, 94, 9596, ■ ■ ■ ) , (65) 

with 

H92^93,9'i,9596,---) = X][n^jJ ^^^'e^ s^^''H92^93,g-i,9596,---) ■ (66) 

We thus see that the wave function has support on holonomies with trivial parallel transport on the cycle, 
92^ 9i^ 9q = 1. There is however one subtlety which is that relating two families of 6j-symbols whose spins 
differ by ±1/2 for example has to be done by hand here. Thus, it may be convenient to restrict to the S0(3) 
model to avoid this. Otherwise, one has to relate the initial conditions of the different families through the 
spin 1/2 recursion [lo| . 

To obtain the above equations, one writes the group Fourier transform of the physical topological state. 



V'phys(5e) = X][n^>] ^Phys(ie) S^^"H9e) , 



{>} 

= E[n^.J '^(Je) E dndnd,. j (-1?^' s^'^Hge) • (67) 

{j'J e' ji,j2,j6 

Using the group averaging identity J dh D^^^l{h)D^^^l{h)D'f^l{h) = l\ (^^^ g ^^^), and the Fourier 

expansion of the Dirac delta on the group, 8{g) = J2j djXjig), the sums over ji,j2,j6 can be explicitly 
performed, 

E ^ndjAe If t f I i-l?^' s^^^H9e) = 6{92' 9i' 9e) s^^'^H92' 93, 9,, 9^96, ■ ■ ■) ■ (68) 

. . . I J4 J5 J6 J 

To get the full projector onto the moduli space of flat connections, one needs the same result on all faces, 
which in general do not have necessarily three links. But it is possible to reduce the analysis to this case by 
performing a sequence of Pachner moves 2-2 to transform a given face with an arbitrary number of links to 
a face with three links. This is feasible since as well-known the 2-2 move is purely kinematical (it just relies 
on SU(2) invariance and not on the dynamics). We shall give a detailed account for this process elsewhere 
since in any case, it would still be necessary to check that the full projector satisfies all the constraints. This 
challenging issue is beyond the scope of the paper. 



V. HIGHER ORDER EQUATIONS 

When using holonomy variables, the quantum theory usually encounters an ambiguity corresponding 
to the choice of the SU(2) representation taken to evaluate these holonomies. Here, to obtain the above 
recurrence relation, we defined the operator H^j using the flux variables as living in the vector represen- 
tation, and the holonomy around the face is naturally in the spin 1 representation. Then, this led to the 
evaluation of the Biedenharn-Elliott identity with a spin being 1. But we could have chosen to express the 
dot products, like X2 ■ Xq in a representation of higher spin, so that the holonomy around the face is in a 
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representation of integral spin, J S N. The question is then how it could affect the dynamical difference 
equation. A natural candidate for the resulting difference equation is the recurrence relation coming from 
restricting the Biedenharn-Elliott identity to spin 2. The latter takes the following form: 

The coefficients look like those of the previous relation, j4_i^o,+1i but for 6j-symbols with a spin 2, 



)^(_ip+i4+i6|i2 i2 21 _^_^2i,+,2+j3+j5+j6('2j, + l) M -^1 H pi (72) 

° ^ ^ ^ ^ be J6 J4j ^ ^ Wi ; J2 Jsj be Je J ^ ^ 

In this section, we describe how to get this relation from the operator Hyj^ by taking it to higher spin. 
This means that we have to perform graspings on spin network nodes in the SU(2) representation of spin 
2. By this, we mean the elementary process which generates the 6j-symbol {-^g } appearing typically in 

(2) 

A};\ As far as we know this technical point has not been discussed in the literature, so we here give some 

details. 

The first natural idea is to square the dot product of the fluxes: {X2 • Xq)'^, but since at the quantum 
level the components X2,X2 do not commute, we take the symmetric part: 

^2 ^2 = 2 (^^2 ^2 + ^2 ^2) = ^2 ^2 + 2 [^2 ' ^2] = ^2 ^2 ~ 2 k ^2 ) (73) 

where in the last step we have used the quantum commutation relation [X^, = —ie'' ^^X!^. We now use 
this symmetrized product to define the dot product of fluxes in the spin 2 representation: 

{X2 ■ Xq)^"^^ = X^^ Xe(j Xgj) - -Xl Xq, (74) 
= {xfXe) ' + I xTXq - ix| Xl . (75) 

The reason for the last term, proportional to is that it is necessary to produce the exact grasping 

for spin 2 as we will show. Indeed, from the previous results (|42|) . (|52|) . we know the action of {X2 ■ ^e)^'^^ 
on a spin network state: 



[X2 ■ X,) sli^^ 



'2^2 \j2j2 l\\ {-iy-+i^+3^ jj2j2 1\ J2(j2 + l)je(je + l) 

^'"^ be h J4i + 2 '"^'^''^^ tie J6 J4. 



(76) 

Parallelly, the Biedenharn-Elliott identity can be used to express the coefficient of the spin 2 grasping, that 
is: l'!-^ f I, in terms of ii^ •'•^ } } and ^ |. The result is nothing but the above factor coming from 
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, , (2) 

the operator (X2 ■ Xq) , so that the latter indeed produces the spin 2 grasping on spin network nodes: 



\Xo • Xr 



(2) 



'tet 



(32 


32 


'11 


lie 


k 





(77) 



Now that we know how to produce 6j-symbols with a spin 2, let us come back to our Hamiltonian operator 
and recurrence relations. In order to generate (j69p . we propose the following operator: 



H, 



(2) 



[Xo-X, 



(2) 



(2) 



6 



(X2-Ad(<72-^5156)(X6)) , 



(78) 



{Xo-X, 



6 



1 



1 



S-(2) Jie} 

^t;,/ «tet 



1 



X2 • Ad(52-l5i56)(X6) + - X2 • - - X2 • Ad(ff2-'<7lff6)(^6). (79) 



(80) 



2 2 

It turns out that this operator gives rise to the expected coefficients, up to an overall factor: 

+2 



fc=-2 



(jl+fcj2,---j6) 
tet 



This can be checked the same way we derive the spin 2 grasping, that is by writing the coefficients A)^ (ji) 
in terms of 6j-symbols with a spin 1, via the Biedenharn-Elliott identity, and compare with the action of 
J. Typically, the Biedenharn-Elliott identity gives: 



ji + 2 ji + 1 1 \ / ji + 1 ii 1 \ _ y^d~Id~^ /ji + 2 ji 2 
32 h is 



i2 i2 is J yjQdj^+iN. 



J2 



32 32 33 



(81) 



This can be lifted to a relation between the coefficients coming from the action of H^ j and those of 



A f ,-,\A ^/dj2-ldj2+ldj^-ldje,+l .(2)..^ 
A+i[ji + l)A+i{]i) = ^TTT^T ^+2\3i)- 



(82) 



Since the left hand side is the coefficient of s. 



(jl+2j2,.--J6) 



(2) 

!>^g^- when acting with H!^ ^ on the spin network function, 

this shows the desired result for the highest shift on ji. The other terms can be derived in the same way. 
The final step is to write the constraint equation: 



E [n(2ie + l)]^(jl,...,j6)4it^=0 



ji,...,j6 e=l 



(83) 



for an arbitrary expansion. By inspection of the coefficients, it comes that: djj_fcyl^],(ji — k) = A*^^ (ji ) , 
so that we can write the following difference equation: 



1(2), 



+2 



E At\3i)^{3i+k,j2, 



,36) 



0. 



(84) 



fc=-2 



However, in contrast with the difference equation coming from H^j, ()63p . the above constraint can not be 
solved with a single initial condition. Indeed, if one looks for instance at the minimum value taken by ji, 

(2) ■ (2) 

that is j™™ = max(|j2 — is|, lis — je\), the coefficients j4_2(if "^), A_^(j™™) become zero. But the equation 
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has three terms, and one has to provide two initial conditions: the values of V'(ji"") aiid ipijf^"^ + 1). So in 
this sense, the spin 2 difference equation is weaker than the spin 1 equation. 

More generally, we can extract from the Biedenharn-Elliott identity a recurrence relation on the 6j-symbol 
of the form: 

for J G N. The coefficients A^^\ji) are products of 6j-symbols with a spin fixed to J, and we expect 
them to admit an expression in terms of 6j-symbols with a spin fixed to 1, through recursive use of the 
Biedenharn-Elliott identity. Typically, the highest coefficient should take the form 

= (-)2^-^+^-^+^3+.5+.a^ jji + J n J I jji + ^1 oc n + k), (86) 

y 32 32 3z\ y JQ 3(i 35) 

up to a normalization depending only on 32, j&- To determine the operator H^'^j producing such an equation, 
one should first find the grasping operator in the representation of spin J. Finally, the corresponding 
Wheeler-DeWitt equation, H^''j V' = 0) will be a difference equation, solved by J initial conditions. 



VI. CONCLUSION 



In this paper, we have introduced a classical constraint H^j on the phase space of loop quantum gravity 
on a single graph. This object naturally comes as a discretization of the scalar constraint H = EfEje'-'^F^f^, 
and is adapted to the specific smearing of the fields at the root of loop quantum gravity. In the holonomy- 
flux variables, it is the projection of the holonomy around a cycle of a spin network graph onto flux variables 
based on a vertex. Prom the point of view of the discrete geometry described in LQG, it amounts to saying 
that the holonomies enable to compute the dihedral angles between triangles as the dihedral angles of flat 
3-simpliccs defined by the lengths. Moreover, 

• it is classically first class (at least on a specific simple triangulation of the sphere), 

• the quantum equation Hyj\il)) = is realized in the spin network basis as a difference equation, 

• on trivalent faces, the latter is a recurrence relation which is known to encode the symmetries at the 
quantum level in the spin foam formalism, 

• and which is equivalent to projecting the face on its flat sector. 

The facts that (in the case we study) the classical algebra is first class and that direct quantization of 
Hyj produces the correct and expected quantum constraints suggest some conclusions about the constraint 
algebra. The usual question of how to split and represent the diffeomorphism algebra in loop quantum 
gravity gets here a simple answer: it is not necessary to split it, but instead to project the curvature along 
all different normals to the triangles meeting at a vertex of the triangulation. Defining the discrete scalar 
constraint iJ^j according to each of these normals is here sufficient. 

To complete the analysis of three-dimensional gravity with those operators, some important things still 
remain to be done. Mainly, one needs to show that it is possible to recover the usual projector onto flat 
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connections for surfaces of arbitrary genus, to find the structure functions of the algebra generated by the 
Hyj classicahy, and to check directly the quantum algebra. Preliminary calculations show that the projector 
is indeed annihilated by all the constraints. But getting an explicit quantum algebra is a difficult task. 

Our result shows that it is possible to precisely identify the recurrence relations satisfied by spin foam 
amplitudes with a quantum implementation of the classical symmetries in loop quantum gravity. This 
suggests two ways to apply this programme to 4d LQG. The first idea is to derive recurrence relations for 
the amplitudes of the most promising spin foam models, and then try to produce them from an operator 



in LQG. A few recurrence relations for non-topological models were obtained in Ul]. Since these models 
describe geometry with areas and normals to triangles, we expect the corresponding operator to produce 
both differences on spins and some differential parts on the normals. 

The second idea is to first derive difference equations from an operator and then interpret them as 
recurrence relations for some spin foam models. Typically, our operator Hyj has a natural generalization in 
4d. There, the curvature tensor has more than one space component on the canonical hypersurface, so we 
may try to sum H^j for three independent faces meeting at a vertex. For example, if is a 4-valent vertex, 
with edges labelled 1, 2, 3, 4, dual to a piece of triangulation, one may consider: 

HvJr2 + Hvj23 + Hv,hi = 0. (87) 

Although certainly too naive, this can be interesting because, if H^ j is seen as the generator of the topological 
symmetry, then the above quantity will only imply invariance under a combination of the corresponding 
changes, as argued in [11] on a toy example. Such an operator would generate a spin foam model, a priori 
different from those recently proposed. 

Ultimately, one may expect these two approaches to coincide. This is actually what we achieved in 
the present article on the 3d model. In four dimensions, some preliminary results have been obtained in 



271 ]. There the topological Ooguri model for BF theory is revisited by lifting the Hamiltonian we have 
just used here, H^j, to the boundary of a 4-simplex. Classically, this provides the phase space of 4d loop 
quantum gravity with a Hamiltonian dynamics which can be interpreted in terms of twisted geometries. In 
the spin network basis, the Wheeler-DeWitt equation gives recursion relations which are actually satisfied 
by the Wigner 15j-symbol. We expect to extend these ideas for more realistic spin foam models than the 
topological model. 



More generally, the relation between spin foams and Hamiltonian dynamics is investigated in 27[ in the 
large spin limit through difference equations on the 4-simplex amplitude. In 3d, such an equation is obviously 
the recursion relation ([T]) on the 6j-symbol, simplified in the asymptotics. Such equations give criteria to 
know whether a model is semi-classically approximated by some quantum Regge calculus. For instance, it 
is clear from these results that the (naive) proposal (j87p is solved in the asymptotics by exponentials of i 
times the Regge action of the 4-simplex with quantized areas. 

Thus, we think that difference and recurrence relations give a very efficient tool to: i) build well-defined 
Hamiltonian operators and ii) make the link with spin foam proposals. It would be interesting to use them 
to make contact and improve recent works which tackle the dynamics of quantum gravity. In particular, a 



new regularization of the Hamiltonian for general relativity has been proposed in 28|]. A relation between 
the topological and some non-topological spin foam models and the canonical framework has been unravelled 
in [2^, but the physical meaning of the projection operators which are used should be clarified. The case 
of spin foams with finite groups has been recently introduced in [30] as a testing ground to investigate such 
issues and also that of the continuum limit. 

That continuum limit issue was not a difficulty here because the theory is topological and the Hamiltonian 
needs not be graph changing, so that the continuum limit could be taken like in [a]. Nevertheless, the point 
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remains challenging for four-dimensional gravity, for which recursion relations on a fixed graph can certainly 
not solve the model. Also in principle, since we have first quantize and then reduce the degrees of freedom 
at the quantum level, the issue of inhomogeneities which plague four-dimensional quantum gravity could be 
addressed in the context of the three-dimensional model, and it is not clear that our quantization provides 
any insight from that view. 

In the framework of group field theory, a symmetry which generates a recursion relation (associated to 
the Pachner move 1-4) in the 3d colored model has been exhibited in 31[. The recursion is a bit different 
from the one we have discussed here (which is in fact associated to the Pachner move 3-2), but the geometric 
interpretation is appealing. Moreover, it could be a way to find generic recursions on spin foam amplitudes 
for larger triangulations than a single simplex. A first step in that direction can be found in [s^], still on a 
simple triangulation of the 3-sphere. 

We want to emphasize that the relation between the Wheeler-DeWitt equation and the Biedenharn- 
Elliott is certainly not a coincidence. Actually the existence of a link between the set of flat connections on 
a surface and Wigner 6j-symbols is known for a while 

BSBQ and more recently p^l] , and our work rather 



sheds light on this intimate relationship. As we have seen in the section HVl the recursion on the 6j-symbol 
and the structure of SU(2) recoupling are strong enough to project the holonomy around a triangular face 
onto the unit of the group. 

Along this line, it may be possible to take advantage of the fact that the recursion relation and the 
BE identity for the homogeneously curved 3d spacetime is known from the quantum group SU(2)q and the 
Turaev-Viro model, in order to find the corresponding operator in LQG. Important results and insights into 
this difficult issue have been recently obtained in 33|], where the non-commutative connection which depends 
on the cosmological constant have been quantized. We think that it could be a good starting point to find 
the Hamiltonian operator which generates the recursion relations for the case of a non-zero cosmological 
constant. 
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